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We show that the interaction of a particle with a directionally solidified interface induces the 
onset of morphological instability provided that the particle-interface distance falls below a critical 
value. This instability occurs at pulling velocities that are below the threshold for the onset of the 
Mullins-Sekerka instability. The expression for the critical distance reveals that this instability is 
manifested only for certain combinations of the physical and processing parameters. Its occurence 
is attributed to the reversal of the thermal gradient in the melt ahead of the interface and behind 
the particle. 
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The freezing of a liquid with a dispersed phase takes 
place in numerous natural and industrial processes. 
Some examples include the formation of ice lenses that 
result from the freezing of soil water Q], the freezing 
of biological cell suspensions in a cryopreservation ex- 
periment the decontamination of metallic pollutants 
from soils 0, the growth of Y123 superconductors by 
the undercooling method and the manufacture of 
particulate reinforced metal matrix composites (PMMC) 
P- The properties of these composite materials are 
enhanced by the addition of the dispersed elements. The 
freezing of a liquid suspension is associated with the 
interaction of the constituents of the dispersed phase 
with a solidifying interface. The first systematic study 
of this interaction was carried out by Uhlmann et al. 
0. They demonstrated the existence of a critical value 
for the growth rate below which the inclusions are 
pushed by the moving interface, and above which they 
are engulfed by the interface and incorporated into the 
solid. Consequently, very low growth rates are conducive 
to particles being pushed by the interface, while high 
growth rates are conducive to particle engulfment. 

The presence of an inclusion in the melt near a 
solid-liquid interface introduces locally a change, albeit 
small, in the thermal gradient ahead of the solid front. 
This, in turn, introduces a small deformation in the 
profile of the interface. The difference in the thermal 
conductivities of the melt and particle stands out as the 
cause for this interfacial deflection 0, H, U • Imagine 
a situation wherein a solid is growing antiparallel to 
the direction of the heat flux and toward an inclusion 
that is less heat conducting than the melt in which it is 
immersed. Then, as the width of the gap separating the 
inclusion from the solid front decreases, heat becomes 
more easily evacuated from the solid phase than from 
ahead, leading to a local reversal of the thermal gradient 
and, consequently, to the local destabilization of the 
interface. As long as the particle remains pushed, 
the disturbance grows and propagates radially with 
decreasing magnitude. 
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FIG. 1: Sketch of a particle of radius a immersed in a melt 
near a deformable solid-liquid interface; aHo is the particle- 
interface distance measured from the particle's center to the 
planar solid front, V is the interface growth rate, hoo is the 
gap width that separates the lowest point on the particle from 
the planar interface, and R represents the radial coordinate 
taken along the planar interface. 



This hand-waving argument is made more precise in 
this paper. The situation considered here is illustrated 
schematically in Fig. 1. An axisymmetric spherical par- 
ticle of radius a is submerged in a binary alloy in such a 
way that the distance from its center to the planar solid- 
liquid interface that emerged from the directional solidi- 
fication of a dilute binary mixture is aHo. We assume a 
microgravity environment. The problem is described by 
the diffusion of solute equation in the melt and the heat 
conduction equations in the melt, solid phase and particle 
in a frame that is moving with velocity V in the upward 
Z— direction; Z is the vertical coordinate (Z > in the 
liquid region) and R is the radial coordinate taken along 
the planar solid-liquid interface. The diffusion of solute 
in the solid phase is neglected. The boundary conditions 
at the solid-liquid interface account for: (i) curvature 
and solute undercoolings in the equation describing the 
equilibrium temperature, (ii) the heat balance, and (iii) 
the conservation of mass. At the particle-melt boundary, 
[Z — (a + /ioo)] 2 + -R 2 = a 2 , we impose the zero mass flux 



2 



and the continuity of the temperature and heat flux; hoo 
denotes the width of the gap between the particle and 
the planar interface. Far away from the solid front, the 
imposed thermal gradients in the solid and liquid phases 
are and G^ s \ respectively, while the concentration 
of solute is maintained at Coo- The nondimensionalisa- 
tion that we adopt makes use of the particle's radius a as 
length scale [Tl|, a/V, the melting point T m , the growth 
rate V , and the concentration far away from the front, 
Cqo, as scales for time, temperature, velocity and con- 
centration, respectively. We have 
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where C is the solute concentration, T is the temper- 
ature, the superscript (q) stands symbolically for par- 
ticle (q — P), solid (q = S) and liquid (q = L), 
= D/DW, where D and D^> are the coefficients 
of solute and thermal diffusion, respectively, e = aV/D 
is the Peclet number, and A r is the Bessel differential op- 
erator, (l/r)d/dr(rd/dr). The corresponding boundary 
conditions are as follows: at the solid-liquid interface, 



T (S) _ T (L) = 1 _ aK + 



Sv n = (kVT (s) - VT {L) ) ■ m 



eC(l - K)v n = -VC • rii 
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At the particle's surface, (z—Hq) 



= 1, the continuity 



of temperature and of the heat flux imply, 

T (i) =r (P) ) ( aVT (P) _ VT< L ') • n P = 0, (6) 

while the zero mass flux condition yields VC ■ np = 0. 
The far field conditions are dT^ L '/dz — » Gl and 
C = 1 as z — > oo, and dT^/dz G s as z — > — oo. 
The symbols that appear in the above equations are 
defined as follows: a = ashjaL is the surface energy 
parameter, where <jsl is the interface excess free energy 
and L is the latent heat of fusion per unit volume; 
M = mCoo/Tm, is the morphological parameter, where 
to is the liquidus slope; S = LVa/(T m k^) is the Stefan 
number; k = k^/k^; a = k^ /k^ L \ where k^ is the 
coefficient of thermal conductance; K is the segregation 
coefficient; v n is the normal growth velocity; and 
np are the unit normal vectors pointing into the melt 
at the interface and particle's surface, respectively; 
G L = aGW/T m and G s = aG^/T m , where is 
the imposed (dimensional) thermal gradient, and k is 
the curvature taken to be positive when the center of 



curvature lies in the soild phase. 

The above set of equations ( with e = 0) 01 admits 
a base state with a planar interface growing at constant 
speed. The determination of the thermal fields in an infi- 
nite medium in which a spherical inclusion has been em- 
bedded has been carried out Here, we make use of 
the method of images 0] to calculate the thermal fields 
in the melt and particle that satisfy the equilibrium tem- 
perature condition at the planar interface, Eq. (3), and 
the conditions at the particle's surface, Eq. (6). Then 
an expression for the thermal field in the solid is derived. 
This expression must satisfy the temperature condition 
at the planar interface, Eq. (3), and must also support 
a planar interface growing with constant speed, Eq. (4). 
We have 



Cb{t, z) = 1, in the melt, 
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T ( B L) {r, z) = l + G L z + M + U(r, z), 
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The plot of the basic profile for the melt's temperature, 
Eq. (8), is shown in Fig. 2. It pertains to an axisym- 
metric particle of radius a = 10/ito whose center is on 
the z— axis and whose separation from the interface is 
0.5/im. It depicts an unstable configuration, wherein 
the thermal gradient in the gap is reversed due to the 
thermally insulating character of the particle. The 
configuration is, however, stable for the case of a highly 
conducting particle. 

In order to examine the stability of the base state, 
we first superimpose axisymmetric, time-dependent in- 
finitesimal disturbances 9^ q \ c and r] upon the basic state 
solutions T B q \ Cb, and the planar interface, respectively. 
We let 



[TM,C\ = [Tg\C B ] + e[eM,c] > 



(12) 



in Eqs. (l)-(6). The resulting equations are then lin- 
earized with respect to the disturbances. The following 
perturbed problem is obtained, 
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A r c+ — = 0, 
oz l 



B {s) = —GsTj + aA r r 



F(r) 



(14) 



6K L ) = -G L r) + aA r r) - F(r)r) + Mc, z = 0, (15) 



rj + Mc, z = 0, (16) 




FIG. 2: Plot of the basic temperature profile in the melt 
(Eq. (8)) for a particle of radius a — 10/im, dimensionless 
gap thickness h x /a — 0.05 and thermal conductance ratios 
a — (left) and a — 10 (right). Note the reversal in the 
temperature gradient in the gap separating the particle from 
the interface for the insulating particle case. The profile is 
evaluated using parameter values for a succinonitrile-acetone 
system (SON- ACE) H0. 
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where the small slope approximation has been invoked in 
the expression for the curvature, k ~ — A r 7y, where rj is 
the perturbation of the planar position, and 
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F(r) = — (r,0) = 2G £ ( — )^ [r2 + (rago)2]6/2 . 



(19) 

The far field conditions for the perturbations are zero 
concentration and zero thermal gradients in the liquid 
and solid phases. The stability problem is solved using 
the finite Hankel transform over the range [0,-^] On 
using the approximation F(r) ~ F(0) as r — > 0, we find 
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where the hat symbol indicates transformed variables, the perturbation f] is obtained. Its solution yields, 
u) is the wavenumber, J m is the Bessel function of the 
first kind of order m and £ stands for the zeta function, 

i-e. C(3)=Er=i"" 3 - A M ' — A( W ) J 



= |M + W(0 ) - iH]e-^/ s 



(23) 



where 



On imposing the heat balance equation at the inter- 
face, Eq. (18), an evolution equation for the growth of 



. / \ cor 2 4C(3)(a-l)Gi 



(24) 
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FIG. 3: The small wavenumber behavior (left) and the large 
wavenumber behavior (right) of the function g(u>) for a sys- 
tem consisting of SCN-ACE containing silicon-carbide (SiC) 
particles. The numerical values of the physical constants are 
[la. 1 16| : liquidus slope m = — 2.8K/wt%, segregation coef- 
ficient K = 1, Coo = 1.3 Wt%, melting point T m = 328 K, 
coefficients of thermal conductivity kS p > = 85Wm -1 K -1 , 
k {s) = 0.225 Wm-'K" 1 and k (L) = 0.223 Wm^K" 1 for 
the particle, solid and liquid, respectively; thermal gradi- 
ents G (S) = 8000 Km" 1 and G (L) = 10800 Km" 1 in the 
solid and liquid, repsectively; a 3 i = 0.009 Jm -2 ; Ho and £ 
have been arbitrarily chosen to equal 1.1 and 10, respectively; 
L = 4.6 x 10 7 Jm~ 3 . 
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Q = kGs + Gl is the conductivity weighted thermal gra- 
dient, and A = + \( L \ The conditions for marginal 
stability are determined by setting the growth rate to 
zero, i.e. A(w) = 0. However, when A = 0, the term 
g{bj)/K{u>) in Eq. (23) is undefined. Therefore, no mean- 
ing to the dispersion relation, A(ui) = 0, could exsit if 
u> did not satisfy the auxiliary equation g(to) = 0. The 
ratio, gf(u>)/A/(u>), where / = djdto is however finite. The 
function g(u>) is singular at uj = and has infinitely many 
zeros Wj, j — 1, 2, 3, . . ., on the postive cu— axis (Fig. 3). 
We set u> = ujj in the dispersion relation A(cj) = and 
solve for Ho . The critical value of the particle-interface 
separation, d c , is the largest value of Hq and is obtained 
by setting u>j equal to the smallest root lo\. We find 



d c . = 



4C(3)(l-o)Gz 



(2 + a)\Q + (1 + k)aul] 
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FIG. 4: (left) Plot of the dimensional critical gap thickness, 
(d c — l)a, as function of the dimensional thermal gradient in 
the melt, G^ L \ in unit of Kelvin per micrometer for selected 
values of the thermal conductance ratio a, (right) plot of the 
critical gap thickness for 2 insulating particles of radii 5fim 
and lO/xm as function of the imposed thermal gradient in the 
meltG (L) . 



Equation (26) describes the instability criterion in terms 
of dimensionless terms. It implies that for prescribed 
thermal gradients in the liquid and solid, and for a given 
thermal conductance ratio a, instability of the planar in- 
terface sets in whenever the dimensionless distance from 
the particle's center to the interface, measured by Hq, 
equals d c . From the physical standpoint, the instabil- 
ity is manifested only for parameter values for which d c 
exceeds unity and the expression (1 — o)Gl is positive. 
Figure 4 depicts the range of parameters for which the 
instability is observable for a system consisting of SCN- 
ACE that is immersed with particles of different sizes and 
thermal conductivities. The plot of the dimensional crit- 
ical gap thickness, = (d c — l)a, as a function of the 
thermal gradient in the melt, G^ L \ for particles of dif- 
ferent thermal conductivities show that the critical gap 
thickness (i) increases with G^ L \ (ii) decreases with a, 
and (iii) particles characterized by high values of a re- 
quire high thermal gradients for the instability to set in. 
For increasing particle's size, the critical gap thickness 
increases as a and G^ remain constant. Furthermore, 
we note that the onset of the instability, being of thermal 
origin, does not depend explicitely on the concentration. 
However, the magnitude of d c is affected by the morpho- 
logical number due to the dependence of the auxiliary 
function g(uj), and thus of the root toi, on M. Depend- 
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ing on the physical process under consideration, the pres- 
ence or absence of this instability can be imposed by the 
selection of an appropriate set of physical parameters as 
dictated by Eq. (26) and shown in Fig. 4. The instability 
put forth in this paper also demonstrates the influence a 
tiny impurity in the melt can have on the morphological 
stability of the interface. 

Ahuja [18j has carried out a detailed experimental study 
of the influence of foreign particles on the morphology 
of a slowly growing solid-liquid interface under normal 



gravity conditions (see also 19]). These experiments de- 
pict in situ the real time evolution of the solid-liquid in- 
terface shape profile in the vicinity of a single particle. 
The sequence of photographs depicting the interaction 
of polystyrene particles in SCN with an initially planar 
solid- liquid interface, shown on page 101 of Ref. [Isj . 
provides qualitative confirmation of the features of the 
instability put forth in this paper. The quantitative val- 
idation of our prediction requires carefully controlled ex- 
periments in a microgravity environment. 
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